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Abstract 

Mesonic fluctuations around the chiral solitons are investigated in the SU(3) chiral quark soliton 
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state. The quantum corrections to the masses (the Casimir energies) are estimated for the octet 
and decuplet baryons. The lack of the confinement in this model demands the cutoff on the energy 
of the fluctuations. Under the assumption that the value of the cutoff energy is 2x(the lightest 
constituent quark mass), these calculation reproduces the masses of the baryons within 15 % error. 
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I. INTRODUCTION 

In the limit of a large number of colors N c Q], QCD reduces to an effective theory of 
the weakly interacting mesons HQ. 

Then the mass of a baryon is proportional to N c 
and the baryons can emerge as solitons of the effective theory [2]. It is widely believed 
that the solitons are important ingredients of the strong interaction at low energies same 
as the chiral symmetry. Some effective models of QCD in low energies region have these 
feature. The Skvrme model |3j is a pure mesonic theory and has_ been investigated in this 
context . The Nambu-Jona-Lasinio (NJL) model 



10, 
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12 and the chiral 



quark soliton model (CQSM) 

,0,0,0,0 

have been studied from a viewpoint of the 
quark model. Usually these models are solved by means of the hedgehog ansatz and the 
cranking method for the non-strange baryons (iV,A). For the hyperons (A,EW,S(*),fi), in 



addition, the embedding of SUY2) meson field to SU(3) one is assumed. 

In the previous works 0,0], we investigated the validity of the hedgehog ansatz for the 
octet and decuplet baryons in the SU(3) CQSM. The isospin vector of the profile function 
takes the non-hedgehog shape for the hyperons and the hedgehog one for the non-strange 
baryon. And the radial component of the profile function for the hyperons approaches to 
the center of the soliton compared with the non-strange baryon. Thus, the shape of the 
soliton changes according to the baryon state. 

These nature of the soliton is due to the quark mass in the body fixed frame. The flavor 
rotation into strange direction rearranges the flavor SU(3) quarks (u,d,s) and couple the 
(u,s) or (d,s) quarks as a SU(2) doublet in this frame. Thus the masses of the doublet are 
asymmetric under a SU(2) transformation. That is reasonable, because SU(2) is the SU(2)v 
or SU(2)u subgroups of SU(3) and the asymmetry is a SU(3) symmetry breaking due to 
the mass difference (m s — m u ) of s and u quarks. The non-hedgehog shape of soliton in our 
approach reflects the asymmetry and the shrink is due to the strange quark mass. 

The approaches using the hedgehog ansatz treat the doublet as SU(2) symmetric com- 
monly to the octet and decuplet baryons and incorporate the SU(3) symmetry breaking by 
a perturbation with respect to the mass difference (m s — m u ). This approach reproduces 
the mass difference between baryons but give the too large absolute values. 

On the other hand, our approach prepares the doublet and the soliton for every baryon 
states and incorporate the SU(3) symmetry breaking at the soliton level in advance. As a 
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result the Hilbert space of the soliton+quark system is enlarged. Our procedure reproduces 



both the absolute values and the mass differences between hyperons [ijj]. However, it gives 
the large absolute values only for the non-strange baryon masses. It is because that there 
is no difference between our approach and other method using the hedgehog ansatz for the 
non-strange baryons. 

For the nucleon mass, there are several attempts to resolve the too large prediction in the 
soliton models. In the Skyrme model, the quantum corrections to the nucleon mass due to 
the meson fluctuations around the hedgehog soliton have been evaluated 120 , 21, 22j. In the 
context of the quark m0 del, the sa m e eo rr eet, ons are e Stlma te d in the NJL mo de, Q Q. 
The corrections give large negative contributions to the non-strange baryon masses, and 
their results are in good agreement with the experimental values. The contributions are 
called the Casimir energies. Their analyses are valid for the non-strange baryons in our 
approach. For the hyperons, however, the non-hedgehog soliton requires some modifications 
on their procedure. 

Thus in this paper we study the Casimir energies due to the meson fluctuations around 
the non-hedgehog soliton and check the consistency of our treatment and the Casimir energy 
for the octet and decuplet baryons. 

In Sec. |n]we show the formal definition of the Casimir energy. In Sec. IIHI after a brief 
review of the SU(3) CQSM, the mean field approximation and the cranking method for 
the non-hedgehog soliton are introduced. Section IIVI is the main part of this paper. At 
first, we define the fluctuations around the soliton and obtain the effective action for the 
fluctuations. Next the normalization and the overlap integral between the fluctuations are 
shown. In Sec. |V]we give the numerical results. Section IVll is the summary of this paper. 

II. CASIMIR ENERGY 

The Casimir energy due to the soliton is given by a difference between the two types of 
the zero point energy: 



AE 

» 3 

where LOi are the energy eigenvalues of the fluctuations around the soliton and denote 
ones in the absence of the soliton. However, this expression diverges quadratically in 3 + 1 
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dimensions. Holzwarth has shown that the ultraviolet divergence requires the three subtrac- 
tiontern^. After the subtractions, the finite expression is given by 



AE 




i^>f |<^)|^(u,f ))f [3 + 6(Vu,f ) 2 - [ Ui /Jp ) 



(2) 



where \z(ui)) and are the energy eigenstates of the fluctuations with the Ui and 



CO 



(0) 



respectively, and (z(uji)\z^ (uo^)} is the overlap integral between these states. This is 
a 3 + 1 dimensional generalization of the result in a 1 + 1 case . 

In this paper we study AE for the octet and decuplet baryons in the CQSM. In the 
Skyrme model j^sj and the NJL model 2s|, AE is dominated by the contributions of the 
zero modes uj~ = and take a negative value for the N and the A. Thus we concentrate on 



the zero modes contributions: 



Wi=0 j 



<^)|# 0) K (0) )> 



(3) 



where X^-=o means that the mode sum for Ui is restricted to the zero modes. 

Equations (J2J) and (J3J) are finite objects for the renormalizable theory in 3 + 1 dimensions. 
In this paper, however we study the CQSM which is a cutoff theory. Thus the sums in 
these equations should be terminated at a physical cutoff point. In addition the overlap 
integrals (z(ui)\z^ (dp)) should be defined by the quantities in this mode. Therefore these 
equations are still formal definitions of the Casimir energy. The cutoff procedure and the 
overlap integrals are defined in Sec. II VI 



III. SU(3) CHIRAL QUARK SOLITON MODEL 

The low energy baryonic states can be obtained through the study of a correlation function 
for the quark operators |2q : 

Mx) = ^e aNc ... a J^ c --- h i) aNjNc {x) . . .i) aih {x), (4) 

J*M = ^ Wo _ A ir-"*V^(y) • ■■r PNjNc (yl (5) 

where N c is the number of color, a and j3 are color indices, both / and g are the spin and 
flavor indices, and is the symmetric tensor representing the baryon state $ in spin and 



4 



flavor space. In the CQSM the correlation function is given by 

1 r Nc 

(Mx)JUy)) = j^jf''"^'^ J DUe^Jlixf.l^-^lyg,), (6) 

where U is the SU(3) chiral meson field and SjrfC/] is the effective action for U. The quark 
Hamiltonian H is given by 

H = -a ■ V + /3 (MW* + rh) , (7) 

where IT 15 = -^-U + ^-^-W, M is the dynamically generated quark mass, and m is the 
current quark mass matrix: 

rh = diag (m u ,m u ,m s ) = m A + m 8 A 8 . (8) 

Then the effective action ^[C/] is defined by the functional determinant for the quark fields, 

iS F [U] = N c log det (id t - H). (9) 
We assume the so-called cranking form 0, for the meson field: 

U^(r,t) = A(t)B\t)W Q 5 (v)B(t)A\t), (10) 

where t/o 5 (r) is the static meson field in the body fixed frame of the soliton, A(t) denotes 
the adiabatic rotation of the system in the SU(3) flavor space, and B(t) describes the spatial 
rotation. Then, the effective action for U 15 is reduced to 

iS F = NAogdet (idt -H'- V A ) , (11) 

where 

V A = -iA^A-iBB\ (12) 

and H' is the rotated quark Hamiltonian with the meson fluctuations defined below. 

For Uq 5 , we assume the embedding of the SU(2) field to the SU(3) matrix but do not 
assume the hedgehog shape jl8|: 

C^»(r) = e n5F( r )A(r) ; ( 13 ) 
where F is the radial component of the profile function, 

A(r) = A ™( r )^ (™ = ±!> 0), (14) 



A m is the contravariant spherical component |28j of a unit vector in isospin space, and A m is 
the covariant one of the SU(2) subalgebra of the Gell-Mann matrices. The transformations 
between the contravariant and the covariant spherical components are given by 

A m = g mm 'A. m , (15) 

m' 

xm = J2g mm '\ m >, (16) 

vn! 

where 

9 mm' = = ( — l) m ^m ,-m' (17) 

We write the flavor rotation 29j as 



■ A(t) , 

.4(/-) = | o y i | .-un. (is) 

where A describes the rotation in SU(2) flavor space and A s represents the rotation into the 
strange directions. Furthermore we parameterize A s (t) as 

V2D(t) \ 

A s (t) = expi\ , (19) 

V2D\t) I 



where D = (Di, D 2 ) T is the isodoublet spinor. In Ref. [3fJ, it is argued that D ~ l/y/~N~ c 
in the large A^ c limit due to the Wess-Zumino term, even if the strange quark mass is light. 
We also employ the classification and treat D perturbatively. 

From the current quark mass matrix (jHJ), we obtain the rotated one: 



rri = A^mA = m X + m 8 D^(A s )X fl , (20) 



where Djj,J(A s ) (//, v = 1, 2, . . . , 8) are the Wigner D functions of A s in the adjoint repre- 
sentation: 

Df}{A s ) = l -tT(A\\,A s \ u ). (21) 

The value of m§ represents the strength of the flavor SU(3) symmetry breaking. 
Here, we define the following quantities: 

k = 2D j D, (22) 
k 3 ee 2D^t 3 D. (23) 



Suppose that the collective variables A and B are quantized and \B) as an eigenstate of a 
collective Hamiltonian. If \B) points to a specific direction in the isospin space, the order 
parameters: 

n m = {B\kq\B), (24) 
k B3 = (B\k 3 \B) (25) 

have nonzero values Then, the expectation value (B\rh'\B) may be approximated by 



rhg — m X + m 8 lim 

tO, 3— »KS0,3 L 



D§(A S )X 3 + D^(A S )X 8 . (26) 



The matrix % is diagonal and its eigenvalues (m Bu ,m B d,m Bs ) are the effective quark 
masses in the body fixed frame. The SU(2) quarks with masses {m Bu ,m B d) interact with 
the chiral field (fT3j). On the other hand, the quark with mg s decouples in the classical 
soliton. 

In Ref. jisl \\\ . we obtain that \kbs\ ~ kbo for the octet and decuplet baryons and the 
value of k bo changes according to the strangeness. For k B3 ps —kbo, the effective quark 
masses become 

m Bu = m u , (27a) 

m B d = rn u cos 2 ^/k^ + m s sin 2 v^so, (27b) 

m Bs = m s + -(to u + m d - 2m s ) sin 2 a/k^. (27c) 
o 



For k B3 ~ kbo. 



m Bu = m u cos 2 v/k^ + m s sin 2 v^so, (28a) 

^Bcz = m u , (28b) 

m Bs = m s + \{m u + m d - 2m s ) sin 2 ^fi^. (28c) 
o 



The important points in the present context are that one of the SU(2) quarks always has 
a light mass m u and the other becomes heavy as k B o grows. In the CQSM, the pseudoscalar 
fields are auxiliary and the corresponding mesons are bound states of the quarks in the body 
fixed frame. Thus, in our treatment, the mesons would be not only the ir but also the K 
and i] according to the strangeness of the soliton. 
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IV. FLUCTUATIONS 



Given Uq 5 , we evaluate the effects of the fluctuations of the pseudoscalar meson fields by 



writing 



21- 



U2 b (r) => Uf (r, t) = e l ^ F{r ^l 2 ■ ■ e ^ F ( r ) A ( r )/ 2 , (29) 

where z(r,t) are the small amplitude fluctuations and have three components in isospin 
space, 

z{r, t) = J2 ^( r ' ( m = ±X > °)' ( 3 °) 

in 

Here, we do not include the genuine SU(3) fluctuations 2^1=4 z<1 ^ a in z(r,t), since we have 
already incorporated the strangeness degrees of freedom by Eq. (jl9j) . The fluctuations z(r, t) 
are induced in the body fixed frame of the soliton and consist of the quarks with the masses 
m Bu and m Bd . Thus, the dispersion relations of the fluctuations also change according to 
the strangeness. 

We expand the effective action up to quadratic order in power of z(r, t). Then the change 
of the chiral field due to the fluctuations is given by 



U? - £/ 75 = E 



E m (r) z m (v, t) - l -P SU 2UT (r) z m (r, t)z m (r, t) 



(31) 



where Psu2 is the projection operator Psu2 = diag (1, 1, 0) for SU(2) subspace and 

E rn (r) = A m (~P S U2 sin F + ry 5 A cos Fj + ry 5 (^\ m - A m A j . (32) 

Using these quantities, the rotated quark Hamiltonian H' with the meson fluctuations in 
Eq. (JTTJ) is given by 

H' = H' + V s , (33) 
H' = ^cx-V + (3 (Mf/ 75 + m B ) , (34) 
V s = (3 (m' - m B ) + (3M {Uf - Uj 5 ) , (35) 

where (rh f — rh B ) determines the fluctuation around the mean field m Bl and H' Q contains the 
effects of the flavor SU(3) symmetry breaking through rh B . 

We expand the correlation function (jUJ) with respect to Vg and Va around the eigenstates 



of the H' Q : 



H' Q \a) = e a \a) (36) 



to obtain the effective action for the collective motions of the soliton and the fluctuations 
around the soliton. In this expansion we use the Schwinger proper time regularization 
13. 123I 31 1 . In this paper, we study only on the fluctuation dependent part of the action. 
The other parts have been studied in Ref. 3] for non-hedgehog soliton. 
The action up to quadratic order of the fluctuations 



dte lu}t z m (r,t), 



(37) 



becomes 



S z = fd 3 r I d 3 r' £ \* m (r, +w) ~z n (r', -u) ( r , r ', u) + g mn 6^ (r - r') (r 

~" " m,n 

I d 3 r^F"(r,^) J m (r,u). (38) 

Here (r) and (r, r', cj) are the local and bilocal kernels respectively: 

(r) = M J2 P R 0«> A ) ( r ) P ^ 2 f/ 75 (r) (r) , (39) 



2^ 



$( 2 ) 



M 2 -°r 2 <*>, A) (r) E m (r) ^ 6 (r) ^ (r) J5„ (r') (r') , (40) 



a.b 



where ipa (r) = (r|a), A is the cutoff parameter, and p K and pr 2 are the cutoff functions shown 
in Appendix |XJ In addition, J m (r,u) are source functions for the fluctuations defined by 



J m (r, w) 



-M ^ p R (e B , A) ^ a (r) £ m (r) </> a (r) 

a 

-Ma" (-w) ^ pr 2 (e a , e 6 , w, A) ^ a (r) £ m (r) V fe (r) (&|/?T», (41) 



0,6 



with 



cr" (u) = J dte^tTXp \rh\t) - m B ] . (42) 
The first term of J m (r, uj) disappears due to the equations of motion for the profile 



function 



14j of the non-hedgehog soliton [l8^ : 

SsinF = (P- A)cosF, 



P m = (P • A)A r 



where 



5(r) = ^p*(e a ,A)4(r)Va(r), 

a 

P m (r) = ^p i? (e a ,A)^ a (r)A m ^ Q (r). 



(43) 
(44) 

(45) 
(46) 



In the second term of J m (r, u>), <r M denote the fluctuations of the mass matrix (|20*|) around 
the mean field states \Bq) of the collective Hamiltonian and approximately satisfy 

(B \& i (-u)\B o ) = 0. (47) 

Thus we finally obtain 

(B o \J m (T,u)\B o ) = 0. (48) 

Up to the same order of approximation as the mean field, the effective action for the 
fluctuations is given by 

Sz= \%\ r dh I dV ^ ¥ m (r ' +u) r (r '' ~ u) ^ (r ' r '' u) + 9mj{3) (r _ r,) (r) ] 

(49) 



m,n 



A. Electric and Magnetic fluctuations 

The kernels $ < - 1 ' ) and $^ 2 - ) are composed of a large number of the quark states and numer- 
ically calculated in our treatment. In this process, the symmetries of the soliton are very 



181 ] is axial symmetric and the profile function 



helpful. The non- hedgehog soliton in Ref. 
satisfies 

F(r,e,V>) = F(r,0,0), (50) 
A m {r,0,ip) = A m (r,9,0)e^. (51) 

Since the third component of the grand spin K(") = J(«) commutes with H' Q [Eq. 
we can select the eigenstate \a) of Hq as the eigenstate of 

K^ ] \a) = M a \a). (52) 

In this case, the isoscalar-scalar S^ a ' h ^ and the isovector-pseudoscalar P„ densities satisfy 

= e W b -M a )v S (a,b) ^ ^ Q) ^ (53) 

P^(r,e,<f) = ^ a (r)i 75 A m ^(r), 

= e i(M b +m-M a )<pp(a,b) ^ ^ q) _ (54) 
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From Eqs. (|5Uj) - (j54j) . the local kernel is independent of the azimuth angle ip, 

$« (r, 9, (p) = (r, 9, 0) , (55) 

and is accurately calculated for the baryon number 5 = 1 case. The part {^) a E m i\)^) of $( 2 ) 
has a simple dependence on (p: 

^ a (r, 0, </?) E m (r, 9, if) (r, 0, ^) = ^)<^ a ( r , 0, ) £ m (r, 9, 0) ^ (r, 5, 0) . (56) 

However it is difficult to calculate $^ 2 - ) itself even numerically, since the bilocal kernel $ ( - 2 * ) 
can contain any pairs of (M a , Mb). Thus, one adopts an ansatz for the fluctuations z and 
simplifies the action (jUJJ). 

In the case of the hedgehog configuration, the fluctuations z are expanded in x m = x m /r 
to separate the radial and angular dependences. According to the parity, the expansion 
branches into two modes : electric z E and magnetic z M fluctuations. Then one usually 



adopt the following ansatz 



3, 



23 



z^ (r, u) = ^2 [xmXrm C Lmi (r,u) + (g mmi - x m x mi ) C Tmi (r, u)] , (57a) 

mi 

z%(r,u) = Yl s mmim2 x m X Mm2 (r,uj). (57b) 

mi ,rri2 

Here, we employ the slightly different parametrization of z^ from Ref. j^J, because of the 
algebraic manageability. 

For the non-hedgehog one, we take the straightforward generalization of Eq. (J57j) . 

z£(t,u) = [AmA mi C Lmi M^)+ (g mmi -A m Ami)C Tmi (^^^)] , (58a) 

mi 

z%(r,u) = Yl e mmim2 k m K Mm2 {r,0^), (58b) 



mi,ni2 
-L,T,M 



where the amplitudes C m depend on r and 9. This ansatz is not so toy as the appearance. 
Expansion of the fluctuations in terms of the vector spherical harmonics is given by 

M = EE \ y jI M 6l(r>»)> (59) 



J,J 3 A=-l 



where (j^(r, co>) is a function of r. By exchanging the order of the summation, we obtain 

z m {r,u) = Xy (J3+m)v5 C/3,m(n 9,u), (60) 
J 3 
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(61) 



where 

J>|J 3 |,A 

Because of Eq. floTj) . Eq. (jSHj) corresponds to Eq. (|5U|) with | J 3 \ < 1. Then the 6* dependences 
of Cm' T ' M correspond to the high J components in Eq. (|6ip. And we can incorporate the 
complicated variations in the direction due to the deformation of the soliton. Thus, Eq. ()58|) 
would work well for low energy excitations, especially the zero modes. 

Using Eq. (J58a|) in Eq. (J49|) and integrating over tp, we obtain the effective action for the 
electric fluctuations: 



Si 



i r du 

2. 2^ 



(2vr) 2 f r 2 sin 9drd9 I r' 2 sin & dr 1 W 

C (r, 61, c) T $%> E (r, 0, r', 5', co) ( Em (r', 9', -co) 

+2tt J r 2 sin 0drd0C™ (r, 0, (r, 0) C* m (r, 0, -w) 



where 



C m M,^) 



C Tm (r, 0, w) 



(62) 



(63) 



and and $( 2 )- E are the 2x2 electric local and bilocal kernels given in Appendix [B] 

respectively. Using Eq. (|58b|) . the effective action for the magnetic fluctuations becomes 



S 



M 



1 I" dco sr^ 

2 2^^ 



r sin 9drd9 / r sin 9'dr'd9 



(64) 



C (r, (9, w) *g>" (r, (9, r', 0', w ) C Mm (r', 0', -a;) 
+2tt y r 2 sin 6drdB(% (r, 5, w) $^ )M (r, 0) ( Mm (r, 9, -co) 

where and are the magnetic local and bilocal kernels given in Appendix E 

respectively. Here the amplitudes ( E ' M with the different spherical indices decouple each 
other. 

The equations of motion for fluctuations (Bethe-Salpeter equations) are obtained from 
the minimization of Sf and Sf' 1 with respect to ( E and £,Jf: 

2tt J r' 2 sin9'dr'd9'^ E (r, 5, r\ 9', u) ( Em (r',9',-u) 

+^ E (r, 9) ( Em (r, 9, -co) = 0, (65a) 
2tt y r' 2 sin 9'dr'd9'^ M (r, 0, r', 6*', w) ( Mm (r', 9', -co) 

+$« M (r, 0) C Mm (r, 0, -co) = 0. (65b) 
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In addition, the boundary conditions 

( Em (R, 9, -u) = 0, (66a) 
( Mm (R,e,-uj) = 0, (66b) 

are required at the surface (r = R) of a sufficient large spherical box, in which the numerical 
calculation is performed. These equations determine the fluctuation amplitudes Cm' M an d 
its eigenvalues ui. 

As referred in Sec. HH we restrict consideration to only zero modes (uij = 0) in B = 1 
sector. Since 

^ E (r, 9, r', 9', 0) = (r', 9', r, 9, 0) , (67) 

(r, 9, r>, 9', 0) = ®<$ M (r', 9', r, 9, 0) , (68) 

we can select the phases of the amplitudes so that 

C™M,0) = C E m (r,9,0), (69) 
&{r,0,0) = ( M m (r,0,0). (70) 

B. Normalization 

The Casimir energy (JHl) needs determining of the normalizations and the overlap integrals 
between the fluctuations in the B = 0, 1 sectors. Weigel et al. define them with the Bethe- 
Salpeter kernels in the NJL model 23j] . Although their calculation is restricted to the case 
of the hedgehog soliton, the procedure can be extended to the non-hedgehog one. For the 
completeness, we repeat their procedure here. At first, we define the metric tensor for each 
spherical index m of the electric and magnetic fluctuations, 

Mt(ry,u Jj ) = £- 2 ^ b (ry i u)l^, (71) 

where &m is the or ' ■ The indices a, b are only for the electric fluctuations and 

refer to the element of 2 x 2 matrix, and ujj are eigenvalues of the fluctuations. Then one 
can demand the normalization condition for the fluctuation amplitudes Cm' M i 

f dW C(r, <*i)Mt (r, r', Uj ) C fem (r', -u s ) = 1. (72) 
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The wave function of the fluctuations are given by 

a ' m (r,^)= /A^( V / Ai)" C (r,x )Wj )r(x, Wj ), (73) 

J c 

where V ' M. m is the square root of the metric tensor determined by 



The normalization condition becomes 

/ d'r^fjr.^rir.w^l. (75) 
Finally the overlap integral between the fluctuations in the B — 0, 1 sectors is defined by 
(UPsWiuf*)) = I rf 3 r^C(r,^)0 ( ° Km (r,^ O) ), (76) 

J a 

where the quantities with the index (0) is in B = sector. 
C. Baryon number B = sector 

In B = sector, it is difficult to calculate accurately the kernels $( 2 ) B ' M ) $( 1 ) £; ' M with the 
spherical plane wave basis 3s| in the large spherical box. The difficulty is remarkable at the 
vicinity of the surface of the box and not negligible since the fluctuations in B = sector 
should freely propagate. 

Instead, we evaluate the dispersion relations of the fluctuations with the help of the 
linear momentum basis. This is systematically done by evaluating the action (J49j) in the 
4-dimensional momentum k space, 

— J2^B q M Bq T(-l,(M Bq /A) 2 ) 
q 



s ' = T / <5fE *•<*)*"<-*) 

M 2 f 1 

-5 m0 —k 2 / d«Vr(0, [M 2 Bq - a(l - a)k 2 ]/A 2 ) 
Stt Jo q 

m 2 r 1 l 

S H i^[k 2 -{m Bu -m Bd ) 2 } J daV (0, [(1 - a)M 2 Bu + aM Bd - a(l - a)fc 2 ]/A^7) 
where q = u,d, and M Bq are the constituent quark masses in the body fixed frame given by 

M Bq = M + m Bq , (78) 
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rriBq are the eigenvalues of the quark mass matrix (|26[). and V is the incomplete gamma 
function. 

By solving the Bethe-Salpeter equation: 

5S* 



0, (79) 



5z m (k) 

we obtain the following dispersion relation, 

k 2 = ( W (0) ) 2 - |k| 2 = Wl (80) 

where W m is the mass of the fluctuation z m (k), uj(°> is the energy eigenvalue, and k is a 
3-dimensional momentum. The mass W m is the lower bound of the lu^ in this sector. The 
Kronecker deltas 5 m o and 5\ m \i in Eq. (fTTjl indicate that the masses of the fluctuations take 
different values according to the spherical index m. An important thing to be emphasized 
in the present context is that the mass W m varies through me, according to the strangeness 
of the soliton. 

On the other hand, the fluctuations are the composite particles and may decay into 
quark- ant iquark pairs at a threshold energy u t h, because of the lack of the confinement 
in this model. If we obey the argument in Ref. [3| as it is and define the threshold so 
that the Feynman parameter (a) integrals in Eq. (fTTj) diverge at k 2 = uf h , we will obtain 
uj t h = 2 min(M Bu , M Bd ) for m = and u th = M Bu + M Bd for m = ±1. In our approach, 
however, there are situations M Bu ^> M Bd or M Bu <C M Bd for the hyperons (Sec. IIII|) . 
Then the "threshold" u t h = M Bu + M Bd for m = ±1 exceeds 2 min(M Bn , M Bd ) largely, 
at which the Dirac sea of the lighter quark may become unstable. Thus the argument in 
Ref. is applicable only for the case M Bu w M Bd . Instead, we employ only one threshold 
ujth = 2 min(M Bu , M Bd ) for m = 0, ±1. Since mm(m Bu , m Bd ) ~ m u for the octet and 
decuplet baryons [Eqs. (|2"7j) . (|2"8j) ]. the fluctuation z m (ci/ )) can be induced within the range: 

W m <uj^ <u th = 2(M + m u ). (81) 

For the wave functions of the fluctuations, we employ the Klein-Gordon operator [2^ as 
the (local +bilocal) kernel in Eq. ()49jl . 

$^(r, r', co) + g mn 5^(r - r')$ (1) (r) = g mn ^\v - v'){u 2 + V 2 - W 2 m ), (82) 

with the ansatz (|57jl . and the boundary condition ()66)1 . Explicit forms of the wave functions 
for B = sector are given in Appendix O The sum in Eq. (jnj) is performed by the 
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TABLE I: The order parameters kbo = |k_b3 I • S is the strangeness of soliton. 

S KBO = \KB3\ 

0.00 

-1 0.80 

-2 1.65 

-3 2.57 



angular momentum (S,P,D-waves) and the momentum (k = |k|) sums. Both S,D-waves are 
the electric fluctuations and P-wave is the magnetic one. The momentum k is discretized 
by the boundary condition (jrJrjj) . The energy region (JHTj) defines the cutoff procedure for the 
zero mode contributions to the Casimir energy. 



V. NUMERICAL RESULTS 



We use the following numerical values for the input parameters: the dynamically gener- 
ated quark mass M = 400 MeV, the current (u,s) quark masses (m„,m s ) = (15,210) MeV, 
and the cutoff parameter A = 700 MeV. These parameters are adjusted to reproduce the 
empirical values of the pion mass and the octet and decuplet baryon masses. 

The order parameters k B o and kbz [Eqs. (J2U) and (J2SJ)] are self-consistently determined 
in the mean field approximation for the octet and decuplet baryons and given in Table UJ 

We first show the amplitudes Cm' T ' M °f the zero mode fluctuations in the B = 1 sector. 
These amplitudes are the solutions for u = of the Bethe-Salpeter (BS) equations (jHHjl and 
are the boundary conditions (JHHJ), and normalized by Eq. (j72)l . Because of the numerical 
reason, we calculate the fluctuation amplitudes multiplied by r. It is enough to evaluate the 
Casimir energy 0. In Fig. [T] the electric = (r(^,r(^) and magnetic r(^f amplitudes 
with the spherical indices m = 0, +1 are displayed for S = 0,-3 cases. The amplitudes 
for S — — 1, —2 have a similar r dependence to one for S = —3. The 9 dependences for 
S — — 1, —2 are somewhat weaker than one for S = —3. 

The amplitudes for S = do not depend on 9. It is consistent with the hedgehog shape 
of the soliton. Since the BS equations and the boundary conditions depend on 9, we can 
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FIG. 1: (Electric and magnetic fluctuation amplitudes) xr with m = 0, +1 in 5 = 1 sector : 
rC4) r Cm' an d r Cm ■ The solid curves represent the 5 = cases, the dashed ones for S = — 3 and 
9 = 0, and the dotted ones for S = —3 and 9 = it/2. 

check the accuracy of the numerical calculation by the spherical symmetry. The curves for 
S 7^ depend not only on r but also 9. The 9 dependence is not direct consequence from 
the non-hedgehog soliton, because we expand the fluctuations z in power of A. However it 
shows the self-consistency of the ansatz (|58|) . 

In our treatment, the meson fluctuations are composed of the quark and antiquark in the 
body fixed frame. Since the quarks for 5* ^ are heavier than ones for S = 0, the fluctuations 
for S 7^ shrink toward the center of the soliton. The same feature applies to the profile 
functions jig ]. The amplitudes, especially , indicate that the meson fluctuations for S ^ 
are quite different from ones for 5 = 0. In the momentum representation, the distributions 
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TABLE II: The masses W m of the meson fluctuations with spherical index m. S is the strangeness 
of soliton. 

S W (MeV) W±i (MeV) 

139 139 

-1 284 311 

-2 315 369 

-3 320 383 

of the fluctuations for S ^ shift to the high energy region comparing with the S = case. 
This nature of the amplitudes is important for the Casimir energy as noted below. 

In the B = sector, the energy eigenvalues of the fluctuations start from the mass 
W m and end at the threshold u t h = 2(M + m u ) = 830 MeV in Eq. (jsTjl . The lower bounds 
W m depend on the spherical index m and the strangeness of the soliton. The values of W m 
are shown in Table |Hj The u quark mass m u is chosen to fit W m with the pion mass for 
5 = 0. It is apparent in this sector that the masses of the fluctuations for S ^ are heavier 
than for S = 0. The masses W + \ and W-\ with the same strangeness are degenerate. It is 
due to the charge (particle-antiparticle) symmetry of the action (|77jl . 

We are now ready to calculate the overlap integrals (z m (u = 0)|zm^(c</ ))) between the 
zero mode fluctuations in B = 1 sector and the S,P,D-waves fluctuations in B = sector. 
These are calculated for individual spherical indices m. Then the radius of the large spherical 
box is chosen as R = 6 fm. Although the overlap integrals and the momenta given by 
Eq. (jMj) depend on the radius, the Casimir energy © do not in principle. In Fig. El we 
present the square of the overlap integrals as a function of It is apparent that the 

D-wave contribution is dominating for electric fluctuations. This result has already been 

nn 

reported for the non-strange baryons in Refs. 22, 23J. 

The sum over in Eq. Q is restricted for S ^ by comparison with S = 0. It is 
because that the overlap integrals for S ^ are shifted to the higher energy region by the 
heavy masses, but the threshold uj t h is in common to the strangeness. 

The Casimir energies AE for S = (0, —1, —2, —3) are summarized in Table IIHI \AE\ 
for S ytz o are about 80 MeV smaller than one for 5 = 0. Thus the heavy fluctuations for 
S are suppressed. It is likely from the physical point of view. However these values are 
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S= 




200 400 600 800 1000 



<o (0) (MeV) 



S = -3 




200 400 600 800 1000 



<d (0) (MeV) 



FIG. 2: The a/ -* dependence of the square of the overlap integral \(z m (u> = 0)\zm ( a; ^)}| 2 with 
m = for the strangeness S = 0, —3. The solid line is the S-wave contribution, the dashed line for 
P-wave, and the dotted line for D-wave. The radius of the large spherical box is chosen as 6 fm. 

affected by the threshold u t h, since the overlap integrals do not sufficiently fall off in the 
vicinity of the uj^ = u t h- In Ref. j^J this phenomenon is found for the electric modes. In 
our approach, we encounter another origin: the heavy fluctuations for the hyperons. For 
example, if we employ u t h = Mb u + Mb s for the fluctuations with the spherical indices 
m = ±1, the Casimir energies for the hyperons become comparable with ones for the non- 
strange baryons. Although such a threshold is rejected from the physical argument in the 
Sec. IIVCI the overlap integrals should exclude the high energy modes normally. However 
it might require the detailed knowledge on the confinement and take us much beyond the 
scope of the present model. Thus if one wants to extract the consistent physical informations 
from the present results, the cutoff procedure for the fluctuations should be regarded as a 
part of the model. 

We summaries the our results in Table HVl Since AE take negative values, the calculated 
masses without AE should be larger than the empirical values. In Ref. 19J, the current s 
quark mass m s = 200 MeV is an adequate value so that the masses without the Casimir 
energies reproduce the empirical ones for the hyperons. Larger m s in this paper pushes up 
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TABLE III: The Casimir energy for B = 1 sector. S is the strangeness of soliton. AE e is the 



contribution of the electric fluctuation. AE rn 


is of the magnetic one. AE = 


A T7> | A 771 

■■ AE e + AE m . 


S AE e (MeV) 




AE m (MeV) 


AE (MeV) 


-73.0 




-222.8 


-296 


-1 -58.2 




-184.7 


-243 


-2 -54.5 




-161.7 


-216 


-3 -53.9 




-160.7 


-215 


TABLE IV: The masses of the octet and decuplet baryons. S is the stranj 


^eness of baryon. E is 


the contribution from the classical soliton and the collective motion. AE 


is tne Casimir energy. 


Expt. is the empirical value. 








S Particle 


E (MeV) 


E + AE (MeV) 


Expt. (MeV) 


N 


1410 


1114 


939 


A 


1648 


1352 


1232 


-1 A 


1224 


981 


1116 


£ 


1249 


1006 






1476 


1233 


1384 


-2 H 


1405 


1189 


1318 




1701 


1485 


1530 


-3 n 


1788 


1573 


1672 



the hyperon masses and cancels out the the Casimir energy. However we can not find the 
moderate values of m s to reproduce the hyperon masses accurately. It is because that the 
soliton breaks for the large m s due to the lack of the confinement in this model. For the 
dynamically generated quark mass M = 400 MeV and the cutoff parameter A = 700 MeV, 
there is the critical value in the vicinity of m s = 210 MeV. As a result, the error due to our 
approach is at most 15% for the octet and decouplet baryon masses. 
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VI. SUMMARY 



We have studied the Casimir energies for the octet and decuplet baryons in the SU(3) 
chiral quark soliton model. They are the quantum corrections due to the meson fluctuations 
around the soliton. Especially we concentrate on the zero mode contributions in the baryon 
number B = 1 sector. Since the soliton takes the non-hedgehog shape for the hyperons 
(A,£(*),H0O,f)) in our approach, we develop a method to treat the fluctuations around the 
non-hedgehog soliton. 

Since the fluctuations are some bound states of the quarks in this model, those equations 
of motion are the Bethe-Salpeter equations. We adopt an ansatz to reduce the degrees of 
freedom of the fluctuations. As a result, the fluctuations are separated into the electric and 
magnetic parts. In addition, each of them can be distinguished with the spherical indices 
due to the axial symmetry of the soliton. Thus we have six independent zero modes. 

The fluctuations for hyperons shrink toward the center of the soliton compared with the 
cases of the non-strange baryons. This feature is because that the quarks constituting the 
soliton become heavy for the hyperons due to the flavor rotation into strange direction. 

The regularized Casimir energy is defined through the overlap integrals between the 
fluctuations in the B = 0, 1 sectors. The integrals can be regarded as the functions of 
the energy eigenvalue of the fluctuations in the B = sector. Then the integrals for the 
hyperons are shifted to the higher energy region than ones for the non-strange baryons. It 
is consistent with the shrink of the fluctuations for hyperons. 

According to the quark picture of the fluctuations, we employ 2x (the lightest constituent 
quark mass) as a cutoff for the Casimir energy. Then the fluctuations for the hyperons are 
suppressed in comparison with the cases of the non-strange baryons. It is reasonable from 
the physical point of view. Our approach reproduces the masses for the octet and decuplet 
baryons within 15 % error at the most. 



APPENDIX A: CUTOFF FUNCTION 



We use the Schwinger proper time regularization 
function for the local density is given by 



11. 



13. 
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in this paper. The cutoff 



p R (e, A) = N cV T + iV c sgn (e) Af H (e, A) 



„val 



(Al) 
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where r)™ 1 is the occupation number of valence quark in the level e and 



sgn(e)^( e ,A) 



dp-^=e~^ = -hgn (e) erfc (|e/A|) 
1/A 2 y/Airp 2 



The cutoff function for the bi-local density is given by 



Pr 2 {e a ,e b ,u}, A) = N t 



uj + e b -e a ' 2 



+ n N cfr 2 {e a ,e b ,u},A) 



where 

/r 2 (e ,e 6 ,a;, A) 
and 



(A2) 



(A3) 



1/A : 



p I e -P e 'i _i_ e -P e 6 , . 

d P\lj^{ " + [uJ 2 -(e a + e b y]R (p,tu,e a ,e b )} , (A4) 



i? (p, e , e fe) = da exp (-p [(1 - e^e 2 , + ae 2 - a(l - a V 2 ] ) . 
Jo 



(A5) 



APPENDIX B: BETHE-SALPETER KERNEL 



By using the ansatz (|58a|) in Eq. ()49)) and integrating over <p, we obtain the action Eq. (J52J) 

/-i \ g 

for the electric fluctuations. The local kernel $m becomes a 2 x 2 diagonal matrix, 



^ E (r, 9) = M 



-l) m A_ m A m 



x 



1 
5* cosF + (P • A) sinF 



-l) m A A 



(Bl) 



where 5 and P m are defined at Eqs. (|43|) . (j4Tj) . Under the situation Eq. (J52J, the quantities 
S 1 , F, (A_ m A m ), and (P ■ A) are independent of <p from the first. 
The bilocal kernel $m is given by 

$V> E (r,e,r',e',u) = (-irM 2 ]Tp r2 (e a ,e^,A) 



a,6 



T 



(a,b) 



Lt a) (r', 6') 
r, B) Lt a) (r', 6K) T<*? 



r,fl)T^' a) (r',0') 
r,#)Ti M (r',0') 



, (B2) 



where 



r>' fc ) (r, 0) = 5 Ma ,M 6+m A m (r, 0, 0) [-5^ (r, 9, 0) sin F + (P^ • A) (r, 0, 0) cos f] (B3) 



t£- 6 > (r, (9) = 5 Ma ,M 6+m ^ L" - A m (r, 5, 0) A" (r, (9, 0)1 (r, (9, 0) 



(B4) 
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and S^'^ and Pm'^ are defined at Eqs. ([53 )1 . ([54 )1 . The integral over (p yields the conservation 
of grand spin 5M a ,M b + m and reduction of computing load. 

For the magnetic fluctuations, by using Eq. (j58b|) . we obtain the action (JMj) . The local 
kernel is given by 

$« M (r,#) = M \l - (-l) m A_ m Aj [,ScosF+ (P- A)sinF] . (B5) 
And the bilocal kernel is given by 

*g>" (r, 9, r', ff, u) = - (-l) m M 2 ^ p V2 (e a , e 6 , W> A) (r, 6») Cg^ (r', 61') , (B6) 

a, 6 

where 



C&> 6 ) (r, 0) = 5 Mo ,m 6+to £ -e m m ^A mi (r, 5, 0) (r, 5, 0) . (B7) 

m\ ,7712 

Similarly, there is the conservation of grand spin SM a ,M b +m- 



APPENDIX C: WAVE FUNCTIONS FOR B = SECTOR 



We define the normalized radial function: 



y/2 

\Ji+i{kR)\ 

where I is the angular momentum, j is the spherical Bessel function, k is the momentum, 
and R is the radius of the large spherical box. The momentum k is determined by the 
boundary condition corresponding to Eq. (fBfij) . 

ji(hR)=0. (C2) 

The normalized wave functions Eq. (|73|) for the electric fluctuations become 

™ I \/(-l) m x^ m x m ( Lm (r,uj) \ 

^ 0)m (r,9,u) = \ Vl ; ^ 1 ; , (C3) 

y y/1- (-l) m x^ m x m ( Tm (r,uj) J 

where ( Lm and ( Tm are the radial part of the wave function, x m = x m /r, and X- m x m is 
independent of (p. For the S-wave fluctuations, the radial parts are given by 

C Lm (r,oo) = -^=u (k r), (C4) 

V47T 



C Tm (r,u;) = -L Uo (fc r), (C5) 

V47T 
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where u = y/W^ + k%. For the D-wave fluctuations, 



C Lm (r,u) = -}=u 2 (k 2 r), (C6) 



( Tm (r-tu) = -J=u 2 (k 2 r), (C7) 

Vo7T 



where a; = y/W^ + k\. 

The normalized wave functions for the magnetic fluctuations become 



^ m (r,9,cu) = y/1 - (-l)™x- m x m { Mm (r,u), (C8) 

where 



t Mm (r,u) = J^u 1 (k 1 r) (C9) 



with = + £3 
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